EXTREMALS FOR THE SOBOLEV INEQUALITY ON THE SEVEN 
DIMENSIONAL QUATERNIONIC HEISENBERG GROUP AND THE 
QUATERNIONIC CONTACT YAMABE PROBLEM 

STEFAN IVANOV, IVAN MINCHEV, AND DIMETER VASSILEV 

Abstract. A complete solution to the quaternionic contact Yamabe problem on the seven dimen- 
sional sphere is given. Extremals for the Sobolev inequality on the seven dimensional Heisenberg 
group are explicitly described and the best constant in the L 2 Folland-Stein embedding theorem 
is determined. 
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1. Introduction 

It is well known that the sphere at infinity of a any non-compact symmetric space M of rank one 
carries a natural Carnot-Carathcodory structure, see [M, P]. A quaternionic contact (qc) structure, 
[Biql, Biq2], appears naturally as the conformal boundary at infinity of the quaternionic hyperbolic 
space. In this paper, following Biquard, a quaternionic contact structure (qc structure) on a real 
(4n+3)-dimcnsional manifold M is a codimension three distribution H locally given as the kernel 
of a R 3 -valued one- form r\ = (771, 772, 773), such that, the three two- forms d,T)i\jj arc the fundamental 
forms of a quaternionic structure on H. This means that there exists a Riemannian metric g on 
H and three local almost complex structures on H satisfying the commutation relations of the 
imaginary quaternions, I1I2I3 = — 1, such that, drjiln = 2<j(/,., .) . The 1-form 77 is determined 
up to a conformal factor and the action of SO (3) on R 3 , and therefore H is equipped with a 
conformal class [g] of Riemannian metrics and a 2-sphere bundle of almost complex structures, the 
quaternionic bundle Q. The 2-sphere bundle of one forms determines uniquely the associated metric 
and a conformal change of the metric is equivalent to a conformal change of the one forms. To every 
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metric in the fixed conformal class one can associate a linear connection preserving the qc structure, 
see [Biql], which we shall call the Biquard connection. 

If the first Pontrijagin class of M vanishes then the 2-sphere bundle of K 3 -valued 1-forms is trivial 
[AK], i.e. there is a globally defined 3-contact form rj that anihilates H , we denote the corresponding 
QC manifold (M, rj). In this case the 2-sphere of associated almost complex structures is also globally 
defined on H . 

Examples of QC manifolds are given in [Biql, Biq2, IMV, Dl]. In particular, any totally umbilic 
hypersurfacc of a quatcrnionic Kahlcr or hyperKahlcr manifold carries such a structure [IMV]. A 
basic example is provided by any 3-Sasakian manifold which can be defined as a (4n + 3)-dimensional 
Riemannian manifold whose Riemannian cone is a hyperKahler manifold. It was shown in [IMV] 
that the torsion endomorphism of the Biquard connection is the obstruction for a given qc-structure 
to be locally 3-Sasakian, up to a multiplication with a constant factor and a S'0(3)-matrix. 

For a fixed metric in the conformal class of metrics on the horizontal space one associates the 
scalar curvature of the associated Biquard connection, called the qc-scalar curvature. Guided by the 
real (Riemannian) and complex (CR) cases, the quaternionic contact Yamabe problem is: given a 
compact QC manifold (M, 77), find a conformal 3-contact form for which the qc-scalar curvature is 
constant. 

In the present paper we provide a solution of this problem on the seven dimensional sphere 
equipped with its natural quaternionic contact structure. The question reduces to the solvability of 
the Yamabe equation (3.4). Taking the conformal factor in the form 77 = u 4: '(Q~ 2 'r), Q = An + 6, 
turns (3.4) into the equation 

Q + 2 2 . , 

Lu = A— — - Au- uScal = — u Seal, 



where A is the horizontal sub-Laplacian, Ah = tr 9 (Vdh), Seal and Seal are the qc-scalar curvatures 
correspondingly of (M, 77) and (M, fj), and 2* = ^p^, with Q = 4n + 6-the homogeneous dimension. 
On a compact quaternionic contact manifold M with a fixed conformal class [77] the Yamabe equation 
characterizes the non-negative extremals of the Yamabe functional defined by 

T(u) = / Ujr^r |Vw| 2 + Scalu 2 W, / u 2 ' dv g = 1, u > 0. 
J M ^ Q — 2 ' Jm 

Considering M equipped with a fixed qc structure, hence, a conformal class [77], the Yamabe constant 
is defined as the infimum 

A(M) = X(M, [77]) = inf{T(u): / u 2 ' dv g = 1, u > 0}. 

Here dv g denotes the Riemannian volume form of the Riemannian metric on M extending in a 
natural way the horizontal metric associated to rj. 

When the Yamabe constant A(M) is less than that of the quaternionic sphere with its standard 
qc structure the existence of solutions to the quaternionic contact Yamabe problem is shown in 
[W], see also [JL1]. We consider the Yamabe problem on the standard unit (An + 3)-dimensional 
quaternionic sphere. The standard 3-Sasaki structure on the sphere is a qc-Einstein structure 77 
having constant qc-scalar curvature Seal = 16n(n + 2). Its images under conformal quaternionic 
contact automorphism have again constant qc-scalar curvature. In [IMV] we conjectured that these 
are the only solutions to the Yamabe problem on the quaternionic sphere. The purpose of this paper 
is to prove this conjecture when the dimension is equal to seven, i.e., n = 1. 

Theorem 1.1. Let fj = j^rj be a conformal deformation of the standard qc-structure fj on the 
quaternionic unit sphere S 7 . If 77 has constant qc-scalar curvature, then up to a multiplicative 
constant rj is obtained from fj by a conformal quaternionic contact automorphism. In particular, 
\(S 7 ) = 48(4tt) 1 / 5 and this minimum value is achieved only by fj and its images under conformal 
quaternionic contact automorphisms. 
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In [IMV] a weaker result was shown, namely the conclusion holds (in all dimensions) provided 
the vertical space of 77 is integrable. We recall the definition of conformal quaternionic contact 
automorphism in Definition 2.1. 

Another motivation for studying the Yamabe equation comes from its connection with the de- 
termination of the norm and extremals in a relevant Sobolev-type embedding on the quaternionic 
Heisenberg group G(H), [GV1] and [Val] and [Va]. As it is well known, the Yamabe equation is 
essentially the Eulcr-Lagrange equation of the extremals for the L 2 case of such embedding results. 
In the considered setting we have the following Theorem due to Folland and Stein [FSt]. 

Theorem 1.2 (Folland and Stein). Let VL a G be an open set in a Carnot group G of homogeneous 
dimension Q and Haar measure dH . For any 1 < p < Q there exists S p = S P (G) > such that for 
u g C °°(fi) 

(1.1) Qf \uf dH(g)j ^ < S p [J \Xu\' dH(g) 

where \Xu\ — Y^JLi \Xju\ 2 with X\, . . . ,X m denoting a basis of the first layer of G and p* = . 

Let S p be the best constant in the Folland-Stein inequality, i.e., the smallest constant for which (1.1) 
holds. The second result of this paper is the following Theorem, which determines the extremals 
and the best constant in Theorem 1.2 when p = 2 for the case of the seven dimensional quaternionic 
Heisenberg group G (H) . As a manifold G (H) = EI x Im H with the group law given by 

(q',U)') = (q„, 0J o ) ° (q,u) = {q a + q,uJ + uj + 2lmq a q), 

where q, q € H and lo,lo € ImH. The standard quaternionic contact(qc) structure is defined by 
the left-invariant quaternionic contact form = (Oi, 2 , 83) = \ (dw — q' ■ dq' + dq' ■ <f ), 
where . denotes the quaternion multiplication. 

Theorem 1.3. Let G (H) = H x LmW be the seven dimensional quaternionic Heisenberg group. 
The best constant in the L 2 Folland-Stein embedding theorem is 

_ 2\/3 
b2 - ^ 

An extremal is given by the function 

v = ^?[(! + M 2 ) 2 + M 2 r 2 > (?,w)eG(H) 

Any other non-negative extremal is obtained from v by translations (5.10) and dilations (5.11). 

Our result confirms the Conjecture made after [GV1, Theorem 1.1]. In [GV1, Theorem 1.6] the 
above Theorem is proved in all dimensions, but with the assumption of partial-symmetry. Here with 
a completely different method from [GV1] we show that the symmetry assumption is superfluous 
in the case of the first quaternionic Heisenberg group. On the other hand, in [IMV] we proved 
Theorem 1.1 in all dimensions, but with the 'extra' assumption of the intcgrability of the vertical 
distribution. In the present paper we remove the 'extra' intcgrability assumption in dimension 
seven. A key step in present result is the establishment of a suitable divergence formula, Theorem 
4.4, see [JL2] for the CR case and [Ob], [LP] for the Riemannian case. With the help of this 
divergence formula we show that the 'new' structure is also qc-Einstein, thus we reduce the Yamabe 
problem on the 7-sphere from solving the non-linear Yamabe equation to a geometrical system of 
differential equations describing the qc-Einstein structures conformal to the standard one. Invoking 
the (quaternionic) Cayley transform, which is a contact conformal diffcomorphism, [IMV], we turn 
the question to the corresponding system on the quaternionic Heisenberg group. On the latter all 
global solutions are explicitly described in [IMV] and this allows us to conclude the proof of our 
results. 
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Remark 1.4. With the left invariant basis of Theorem 1.3 the Heisenberg group G (H) is not a 
group of Heisenberg type. If we consider G (H) as a group of Heisenberg type then the best constant 
in the L 2 Folland-Stein embedding theorem is, cf. [GV1, Theorem 1.6], 



tt 2 /5 2V2' 

and an extremal is given by the function 

F(q,u) = 7 [(1 + M 2 ) 2 + IBM 2 )]" 2 , (P)£G(H) 

where 

7 = 32^- 17 / 50 2 1 / 5 15 2 / 5 . 

Organization of the paper. The paper uses some results from [IMV]. In order to make the 
present paper self-contained, in Section 2 we give a review of the notion of a quaternionic contact 
structure and collect formulas and results from [IMV] that will be used in the subsequent sections. 

Section 3 and 4 are of technical nature. In the former we find some transformations formulas for 
relevant tensors, while in the latter we prove certain divergence formulas. The key result is Theorem 
4.4, with the help of which in the last Section we prove the main Theorems. 

Convention 1.5. We use the following conventions: 

• {ei, . . . ,e4„} denotes an orthonormal basis of the horizontal space H. 

• The summation convention over repeated vectors from the basis {ei, . . . , e4„} will be used. 
For example, for a (0,4)-tensor P, the formula k = P(eb, e a , e a , &b) means 

in 

k = ^ P(e b ,e a ,e a ,e b ). 

a, 6=1 

• The triple (i,j,k) denotes any cyclic permutation of (1,2,3). 

Acknowledgements S.Ivanov is visiting Max-Plank-Institut fur Mathematics, Bonn. S.I. thanks 
MPIM, Bonn for providing the support and an excellent research environment during the final stages 
of the paper. S.I. is a Senior Associate to the Abdus Salam ICTP. I.Minchev is a member of the 
Junior Research Group "Special Geometries in Mathematical Physics" founded by the Volkswagen 
Foundation. The authors would like to thank The National Academies for the financial support and 
University of California, Riverside and University of Sofia for hosting the respective visits of the 
authors. 

The authors would like to thank the referee for remarks making the exposition clearer and 
spotting several typos in the paper. 

2. Quaternionic contact manifolds 

In this section we will briefly review the basic notions of quaternionic contact geometry and recall 
some results from [Biql] and [IMV]. 

For the purposes of this paper, a quaternionic contact (QC) manifold (M, g, Q) is a An + 3 
dimensional manifold M with a codimension three distribution H equipped with a metric g and an 
Sp(n)Sp(l) structure, i.e., we have 

i) a 2-sphcrc bundle Q over M of almost complex structures, such that, we have Q = {al\ +6/2+ 
c/3 : a 2 + b 2 + c 2 = 1}, where the almost complex structures I s : H — > H, I 2 = — 1, s = 
1, 2, 3, satisfy the commutation relations of the imaginary quaternions — —I2I1 = I3; 

ii) H is the kernel of a 1-form r\ = (771 , 772, 773) with values in K 3 and the following compatibility 
condition holds 



2g(I s X,Y) = dn s (X,Y) 7 a = 1,2, 3, X, Y e H. 



EXTREMALS FOR THE SOBOLEV INEQUALITY 



5 



Correspondingly, given a quaternionic contact manifold we shall denote with 77 any associated 
contact form. The associated contact form is determined up to an 50(3)-action, namely if ^ e 5*0(3) 
with smooth functions as entries then tyrj is again a contact form satisfying the above compatibility 
condition (rotating also the almost complex structures) . On the other hand, if we consider the 
conformal class [g] , the associated contact forms are determined up to a multiplication with a positive 
function \x and an 50(3)-action, namely if W G 50(3) then n^fr] is a contact form associated with 
a metric in the conformal class [g]. 

We shall denote with (M, 77) a QC manifold with a fixed globally defined contact form. A special 
phenomena here, noted in [Biql], is that the 3-contact form 77 determines the quaternionic structure 
and the metric on the horizontal bundle in a unique way. 

A QC manifold (M,g,<Q) is called conformal to (M, <?,Q) if g G [g]. In that case, if 77 is a 
corresponding associated one- form with complex structures I s , s = 1,2,3, we have 77 = /j<J/77 for 
some $ G 50(3) with smooth functions as entries and a positive function fi. In particular, starting 
with a QC manifold (M, 77) and defining 77 = fiij we obtain a QC manifold (M, 77) conformal to the 
original one. 

Definition 2.1. A diffeomorphism 4> of a QC manifold {M, [g],Q) is called a conformal quaternionic 
contact automorphism (conformal qc-automorphism) if <fi preserves the QC structure, i.e. 

cj)*r] = /Lt^7 • 77, 

for some positive smooth function [i and some matrix W G 50(3) with smooth functions as entries 
and rj = (771, 772, 773)* is a local 1-form considered as a column vector of three one forms as entries. 

Any endomorphism 'J of H can be decomposed with respect to the quaternionic structure (Q, g) 
uniquely into 5j)(n)-invariant parts as follows = >]/+++ + \]H + v|> 1 + v|> h , where 
\I>+++ commutes with all three h, commutes with I\ and anti-commutes with the other two 

and etc. The two 5p(7i)5p(l)-invariant components are given by 

(2.1) * [3] =*+++, *[_!]= * + — + *-+- + 

Denoting the corresponding (0,2) tensor via g by the same letter one sees that the Sp(n)Sp(l)- 
invariant components are the projections on the eigenspaces of the Casimir operator 

(2.2) f = h ® h + h®h + h®h 

corresponding, respectively, to the eigenvalues 3 and —1, see [CSal]. If n = 1 then the space of 
symmetric endomorphisms commuting with all Ii, i = 1, 2, 3 is 1-dimensional, i.e. the [3]-component 
of any symmetric endomorphism ^ on H is proportional to the identity, ^[ 3 ] = tr Idijj. 

On a quaternionic contact manifold there exists a canonical connection defined in [Biql] when 
the dimension (4n + 3) > 7, and in [D] in the 7-dimcnsional case. 

Theorem 2.2. [Biql] Let (M,g,Q) be a quaternionic contact manifold of dimension 4n + 3 > 7 
and a fixed metric g on H in the conformal class [g] . Then there exists a unique connection V with 
torsion T on M 4 "+ 3 and a unique supplementary subspace V to H in TM , such that: 

i) V preserves the decomposition H © V and the metric g; 
ii) for X,Y G H, one has T{X,Y) = -[X,Y]\ V ; 

Hi) V preserves the Sp(n)Sp (1) -structure on H , i.e., V<? = and VQ C Q; 
iv) for ^ G V, the endomorphism T(£, .)| ff of H lies in (sp(n) © sp(l)) _L C gl(4n); 
v) the connection on V is induced by the natural identification (p of V with the subspace sp(l) of 
the endomorphisms of H , i.e. V</2 = 0. 

We shall call the above connection the Biquard connection. Biquard [Biql] also described the 
supplementary subspace V explicitly, namely, locally V is generated by vector fields {Cij^jCa}) 
such that 

^ 23 ^ Vs((,k) = Ssk, {£ s ->dr) s )\ H = 0, 

(£ s jd77 fe )| H = -(£k^dn s )\ H . 
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The vector fields £1,^2, £3 are called Reeb vector fields or fundamental vector fields. 

If the dimension of M is seven, the conditions (2.3) do not always hold. Duchcmin shows in 
[D] that if we assume, in addition, the existence of Reeb vector fields as in (2.3), then Theorem 2.2 
holds. Henceforth, by a qc structure in dimension 7 we shall mean a qc structure satisfying (2.3). 

Notice that equations (2.3) are invariant under the natural SO(3) action. Using the triple of Reeb 
vector fields we extend g to a metric on M by requiring span{^\, £2, £3} = V _L H and g(£ s , £&) = 5 s k- 
The extended metric docs not depend on the action of 50(3) on V, but it changes in an obvious 
manner if r\ is multiplied by a conformal factor. Clearly, the Biquard connection preserves the 
extended metric on TM, Vg = 0. We shall also extend the quaternionic structure by setting I s \ v = 0. 
The fundamental 2-forms cjj, i = 1,2,3 of the quaternionic structure Q are defined by 

(2.4) 2u ilH = di lAH , £jwi = 0, ieV. 

Due to (2.4), the torsion restricted to H has the form 

3 

(2.5) T(X,Y)=-[X,Y]\ v =2'£w a (X,Y)Z t , X,Y e H. 

s=l 

The properties of the Biquard connection are encoded in the properties of the torsion endomor- 
phism T% = T(£, .) : H — > H, £ e V. Decomposing the endomorphism T% G (sp(n) + sp(l)) 1 - into 
its symmetric part T® and skew-symmetric part b^, = T® + b^, we summarize the description of 
the torsion due to O. Biquard in the following Proposition. 

Proposition 2.3. [Biql] The torsion is completely trace-free, 

trTt =s(T 5 (e„),e„) = 0, trT c o I = g (T ( (e a ),Ie a ) = 0, I € Q, 

where e\...e± n is an orthonormal basis of H . Decomposing the torsion into symmetric and anti- 
symmetric parts, = Tj} ( + i = 1,2,3, we have: the symmetric part of the torsion has the 
properties 

Tlh = -hTl 

h{Tl) + - = h{Tl)- + - , h{Tl)-+- =I 2 {Tl)-+, h{Tl)-+ =h{Tl)+-; 
the skew- symmetric part can be represented in the following way 



where u is a traceless symmetric (l,l)-tensor on H which commutes with I\, 12,13 ■ 

If n = 1 then the tensor u vanishes identically, u = and the torsion is a symmetric tensor, 
T c = T°. 

The covariant derivative of the quaternionic contact structure with respect to the Biquard con- 
nection and the covariant derivative of the distribution V are given by 

(2.6) VJi = -a j ® I k +a k ® Ij , V& = -a, ® & + & k ® & , 
where the sp(l)-connection 1-forms a s on H are given by [Biql] 

(2.7) a i (X) = dr) k (£ j ,X) = -di lj (Z k ,X), X e H, & e V, 

while the sp(l)-connection 1-forms a s on the vertical space V are calculated in [IMV] 

f Seal 1 

(2.8) = <Mfr,&)- S is ( 16n(n + 2) + 2 + + <%(&,6)) 

where s € {1, 2, 3}. The vanishing of the sp(l)-connection 1-forms on H is equivalent to the vanishing 
of the torsion endomorphism of the Biquard connection, see [IMV] . 
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2.1. The qc-Einstein condition and Bianchi identities. We explain briefly the consequences 
of the Bianchi identities and the notion of qc-Einstein manifold introduced in [IMV] since it plays a 
crucial role in solving the Yamabe equation in the quaternionic seven dimensional sphere. For more 
details see [IMV]. 

Let R = [V, V] — V[ , ] be the curvature tensor of V. The Ricci tensor and the scalar curvature 
Seal of the Biquard connection, called qc-Ricci tensor and qc-scalar curvature, respectively, are 
defined by 

Ric(X, Y) = g(R(e a ,X)Y, e a ), X,Y £ H, Seal = Ric(e a , e a ) = g{R{e b , e a )e a , e b ). 

According to [Biql] the Ricci tensor restricted to H is a symmetric tensor. If the trace- free part of the 
qc-Ricci tensor is zero we call the quaternionic structure a qc-Einstein manifold [IMV]. It is shown 
in [IMV] that the qc-Ricci tensor is completely determined by the components of the torsion. First, 
recall the notion of the iSp(n)Sp(l)-invariant trace- free symmetric 2-tensors T°, U on H introduced 
in [IMV] by 

T°(X, Y) d = f g({Tlh + T° 7 2 + T° J 3 )X, Y), U(X, Y) = f g(uX, Y), X,Y e H. 

The tensor T° belongs to [— l]-eigenspace while U is in the [3]-eigenspace of the operator f given by 
(2.2), i.e., they have the properties: 

(2.9) T°(X, Y) + T a (hX, hY) + T°(I 2 X, I 2 Y) + T°(I 3 X, I 3 Y) = 0, 

(2.10) 3U(X, Y) - U(hX, hY) - U(I 2 X, I 2 Y) - U{I 3 X, I 3 Y) = 0. 

Theorem 1.3, Theorem 3.12 and Corollary 3.14 in [IMV] imply: 

Theorem 2.4. [IMV] Let (M 4n+3 , g, Q) be a quaternionic contact {An + 3)- dimensional manifold, 
n > 1 . For any X,Y £ H the qc-Ricci tensor and the qc-scalar curvature satisfy 

Ric(X,Y) = (2n + 2)T°{X,Y) + (4n+lO)U(X,Y) + ^g(X,Y) 
Seal = -8n(n + 2) 5 (T(£i,6U 3 ) 

For n = 1 the above formulas hold with U = 0. 

In particular, the qc-Einstein condition is equivalent to the vanishing of the torsion endomorphism 
of the Biquard connection. If Seal ^ the latter holds exactly when the qc-structure is 3-Sasakian 
up to a multiplication by a constant and an SO(3)-matrix with smooth entries. 

For the last part of the above Theorem, we remind that a (4n+3)-dimcnsional Ricmannian man- 
ifold (M,g) is called 3-Sasakian if the cone metric gw = t 2 g + dt 2 on N = M x R+ is a hyperkahlcr 
metric, namely, it has holonomy contained in Sp(n + 1). 

The Ricci 2-forms p s , s = 1, 2, 3 of a quaternionic contact structure are defined by 

An p s (B, C) = g(R(B, C)e a , I s e a ), B,C £ T(TM). 

For ease of reference, in the following Theorem we summarize the properties of the Ricci 2-forms, 
the scalar curvature and the torsion evaluated on the vertical space established in Lemma 3.11, 
Corollary 3.14 Proposition 4.3 and Proposition 4.4 of [IMV] . 

Theorem 2.5. [IMV] The Ricci 2-forms satisfy 
P i(X,Y) = 2g{{Tl)—+hX,Y)- 

(2.11) P2 (X,Y) = 2g((Tl) + —hX,Y)- 
p 3 (X,Y) = 2g{{Tl)- + -I 2 X,Y)- 



-2g(huX, Y) 

- 2g{I 2 uX, Y) 

- 2g(I 3 uX, Y) 



Seal 



8n(n + 2, 

Seal 
8n(n + 2) 

Seal 
8n(n + 2) 



to 3 (X,Y). 
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(2.12) Pi(X,Z j )=w j ([Z J ,ik],X), Pi (X,a)=co k ([^,^},X), 
Pi(I k X,^) = -fH(I j X,S k )=g{T(t j ,Z k ),I i X)=u> i ([Z j ,Z k ],X), 

(2.13) Pi(Zi,Zi)+Pk(Sk,Zj) = 8n fn + 2) ^ Scal) - 

The torsion of the Biquard connection restricted to V satisfies the equality 

(2.14) rte>6) = __^_ & _ [6) ^ i 

where [£i,£j]H denotes the projection on H parallel to the vertical space V. 

We also recall the definition of the Sp(n)Sp(l)-uxvaxiaxit vector field A, which appeared naturally 
in the Bianchi identities investigated in [IMV] 

A = Ii[fc,6]+/aK3 J fc]+/3Ki,&]. 

We shall denote with the same letter the corresponding horizontal one-form, i.e., 

A(X) = g(Ii[b,&]+h[k,Zi] + h[Zi,b],X). 

The horizontal divergence V*P of a (0,2)-tensor field P on M with respect to Biquard connection 
is defined to be the (0,l)-tensor field 

V*P(.) = (V eo P)(e ,.). 
Then we conclude from [IMV, Theorem 4.8], that 

Theorem 2.6. [IMV] On a {An + 3)- dimensional QC manifold with constant qc-scalar curvature we 
have the formulas 

i _ , 7 

(2.15) V*T° = (n + 2)A, V*U = — - — A 

3. CONFORMAL TRANSFORMATIONS 

Note that a conformal quaternionic contact transformation between two quaternionic contact 
manifold is a diffeomorphism $ which satisfies 

$*7y = n VP • t], 

for some positive smooth function /i and some matrix \& S 50(3) with smooth functions as entries 
and r\ is an R 3 -valued one form, r\ = r/2, 7/3)' is a column vector with entries one-forms. The 
Biquard connection does not change under rotations, i.e., the Biquard connection of \& • 77 and 77 
coincide. Hence, studying conformal transformations we may consider only transformations &*r/ = 
/j 77. 

Let ft, be a positive smooth function on a QC manifold (M,n). Let fj = 7^77 be a conformal 
deformation of the QC structure 77. We will denote the objects related to fj by over-lining the same 
object corresponding to rj. Thus, dfj = — dhArj + j^dr] and g = ^-<7. The new triple {£i,£2,£3} 
is determined by the conditions defining the Reeb vector fields. We have 

(3.1) & = 2h£ s + I s Vh, s= 1,2,3, 

where V/i is the horizontal gradient defined by g(Vh, X) = dh(X), X E H. 

The components of the torsion tensor transform according to the following formulas from [IMV, 
Section 5] 

(3.2) T°(V,V) = T°(X,Y) + h-^Vdh^^X^), 

(3.3) U{X,Y) = U(X,Y) + (2h)- 1 [Vdh-2h- 1 dh®dh} [3][0] (X,Y), 
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where the symmetric part is given by, cf. (3.9), 

3 

[Vdh] [sym] (X,Y) = Wdh(X,Y) + J2 dh (Zs)us(X,Y) 

s=l 

and [3][o] indicates the trace free part of the [3]-component of the corresponding tensor. In addition, 
the qc-scalar curvature changes according to the formula [Biql] 

(3.4) Seal = 2h (Seal) - 8(n + 2) 2 h^\Vh\ 2 + 8(n + 2) Ah. 
The following vectors will be important for our considerations, 

(3.5) Ai = ii hence A = A x + A 2 + A 3 . 

Lemma 3.1. Let h be a positive smooth function on a QC manifold (M,g,<Q) with constant qc- 
scalar curvature Seal = 16n(n + 2) and fj = ^- r\ a conformal deformation of the qc structure rj. If 
fj is a 3-Sasakian structure, then we have the formulas 

(3.6) Ax{X) = ~h~ 2 dh(X) - hi^\Vh\ 2 dh{X) 

- h-^wdhihx,^) + wdh(i 3 x,^+ hi- 2 (dh{&)dh(i 2 x) + dh(&)dh(i 3 x) 

+ ^h- 2 ^Vdh{I 2 X,I 2 ^h) + Vdh{I 3 X,I 3 Vh) 

The expressions for A 2 and A 3 can be obtained from the above formula by a cyclic permutation of 
(1,2,3). Thus, we have also 

A{x) = -hi- 2 dh{x) - hi-3\vh\ 2 dh{x) 



3 3 -.3 

h- 1 ^ Vdh(I s X,£ s ) + h- 2 Y,dh{i s )dh{I s X) + ^h- 2 Y, Vdh{I s X,I s Vh) 



s=l 



Proof. First we calculate the sp (l)-connection 1-forms of the Biquard connection V. For a 3- 
Sasaki structure we have dfji(t;j,t;k) = 2, £;jg??7« — 0, the non-zero sp(l)-connection 1-forms are 
<3i(£i) = — 2,i = 1,2,3, and the qc-scalar curvature Seal = 16n(n + 2) (see [Example 4.12,[IMV]]). 
Then (3.1), (2.7), and (2.8) yield 

2<%fe,&) = 2II- 1 + h~ 2 \\dhf, on(X) = -h- x dh{hX), 

Oifc) = -h^dhfa) = -ctjfa), 4ai(6) = -4 - 2/j- 1 - h- 2 \\dh\\ 2 . 

From the 3-Sasakian assumption the commutators are = — 2£&. Thus, for X G H taking 

also into account (3.1) we have 

g([Zi,b],hX) = -2gfa,hX) = -2g(2h& + I 3 Vh,I 3 X) = -2dh(X). 

Therefore, using again (3.1), we obtain 

(3.8) -2dh(X) = g([£i, IaX) = g([2h^ + /iVh,2hf a + I 2 Vh], I 3 X) 
= -4h 2 A 3 (X) + 2hg{[{ 1 , I 2 Vh], I 3 X) + 2hg{ [I x Vh, &], I 3 X) 

+ g([hVh, I 2 Vh],I 3 X). 

The last three terms are transformed as follows. The first equals 

g([^, I 2 Vh],I 3 X) = ff((V Cl J 2 )Vft + / 2 V ?1 V/i, I 3 X) - g I 2 Vh), I 3 X) 

= -a 3 (b)dh{I 2 X) + ai (^)dh(X) - VdhituhX) - g(T(^,I 2 Vh),I 3 X), 
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where we use (2.6) and the fact that V preserves the splitting H © V. The second term is 

g{[hVh,i 2 lhX) = a 2 (&)dh(X) + a 3 (&)dh(hX) - Vdh(&,I 2 X) 

- g(T(hVh,b),I 3 X), 

and finally 

g{[hVh,hVh},I 3 X) = -a 3 (hVh)dh(I 2 X) + ai (hV h) dh(X) - Vdh(hVh, hX) 

+ a 2 {I 2 Vh)dh{X) + a 3 (I 2 Vh)dh(hX) - Vdh(I 2 Vh,I 2 X). 

Next we apply (3.7) to the last three equalities, then substitute their sum into (3.8), after which we 
use the commutation relations 

3 

Vdh(X,Y)-Vdh(Y,X) = -dh(T(X,Y))= -2 5> a (X,Y)dfc(6). 

Vdh(X,£) -Vdh(£,X) = -dh(T(X,£)), X, Y e H, S,eV. 
The result is the following identity 

(3.10) Ah 2 A 3 {X) = (-4/i + h^WVhW^dhiX) 

-2h[Vdh[I 1 X,£ 1 ) + \7dh{I 2 X,&)]- [Vdh (I X X, hVh)) + Vdh (I 2 X, I 2 Vh)] 
+ 2[dh(£ 1 )dh(I 1 X) + dh(&) dh(I 2 X) + 2 dh(&) dh(I 3 X)] 
+ 2h[Tfa,hX,Vh) + T(&,I 2 X,Vh) - T(h,I 2 X,I 3 Vh) + T(£ a , hX, I 3 Vh)], 

where T(£, X, Y) = g(T^X, Y) for a vertical vector £ and horizontal vectors X and Y. With the help 
of Proposition 2.3 we decompose the torsions into symmetric and anti-symmetric part 2^ = T°. +IJJ ', 
i = l, 2, 3, and then express the symmetric parts of the torsion terms in the form 2^ = (2^) + 
(T°)-+-, T° - (T°)— + + (T°)+~ . Hence, using T°~ + = 2(2*)+— 1 2 = 2(1*)"+-/! etc., 
which follows again from Proposition 2.3, the sum of the torsion terms in (3.10) can be seen to equal 
2T°~~ + (X, V/i) - AU{X,Vh). This allows us to rewrite (3.10) in the form 

(3.11) AA 3 (X) = {-Ah' 1 + h- 3 \\Wh\\ 2 )dh(X)- 2 hT 1 [ Vdh (I X X, &) + Vdh(I 2 X,&)] 

+ 2h- 2 [dh(£ 1 )dh{I 1 X) + dh(&)dh(I 2 X) + 2dh(&)dh(I 3 X)] 
- hT 2 [Vdh(hX,IxVh) + Vdh(I 2 X,I 2 Vh)] + Ah,- 1 [ (T°" + (V/i, X) - 2U(Vh,X)]. 
Using (3.2) the T° + component of the torsion can be expressed by h as follows, see (2.1) and (2.9), 

4T°~ + (Vh,X) = T°(Vh,X) ~ T°{hVh,hX) - T°(I 2 Vh,I 2 X) + T°(I 3 Vh,I 3 X) 
= -/i- 1 {[Vd/i] hl] (V/i,X)-[V^] hl] (/iV/ l ,/iX)-[V^] [ _ 1] (/ 2 V/ l ,/ 2 X) + [V^] [ _ 1] (/ 3 V/i,/3^)} 

3 

-hT 1 53{dfc(£») [g(I.Vh,X) - 9 {I s hVh,hX)- g(I s I 2 Vh,I 2 X) + g(I a I a Vh,hX)]} 

s=l 

= -h- 1 (vdh{Vh,X) - Vdh(I{Vh,IiX) - Vdh(I 2 Vh,I 2 X) + Vdh(I 3 Vh,I 3 X)\ 

+ Ah- 1 dh(&) dh(I 3 X). 

Invoking equation (3.9) we can put V/i in second place in the Hessian terms, thus, proving the 
formula 

(3.12) 4T°~ + (V/i,X) = - Ah' 1 dhfe) dh(I 3 X) 

- h- 1 tvdh(X,Vh) - Vdh(hX,hVh) - V dh (I 2 X , I 2 V h) + Vdh (I 3 X, I 3 Vh) }. 
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On the other hand, (2.10), (3.3) and the Yamabe equation (3.4) give 

3 

(3.13) 8U(Vh,X) = -h- 1 {vdh(Vh,X) + Vdh(I a Vh,I s X) 

s=l 

-2h- 1 \\Vh\\ 2 dh(X) - —dh(X) + 2h~ 1 ^^- dh(X)\ 

n n J 

3 

= -h' 1 {vdh{Vh,X) + ^2Vdh(I s Vh,I s X)\ 



8 = 1 

,-1 



- h-*l-2h-Wdh(X) - *»-4nh + {n + 2)h-i\\VhV + ^ jW | 

I n n J 

3 

= -/i" 1 |vd/i(X,V/i) + ^Vcto(I s X,I s V/i)} - /i -1 (-3/i -1 ||V/i|| 2 - 2 + Ah)dh{X). 

8=1 

Substituting the last two formulas in (3.11) gives A3 in the form of (3.6) written for A\, cf. the 
paragraph after (3.6). □ 

4. Divergence formulas 

We shall need the divergences of various vector/forms through the almost complex structures, 
so we start with a general formula valid for any horizontal vector/form A. Let {e\, . . . , e4„} be an 
orthonormal basis of H . The divergence of I\A is 

V*(hA) = (V ea (hA))(e a ) = -(V ea A)(he a ) - A{{V e J x )e a ), 

recalling hA(X) = -A(hX). 

We say that an orthonormal frame 

{ei,e 2 = /iei,e 3 = I 2 ei,e4 = hex, . . . , e 4 „ = hdnS, £1, 6, £3} 
is a qc-normal frame (at a point) if the connection 1-forms of the Biquard connection vanish (at that 
point). Lemma 4.5 in [LMV] asserts that a qc-normal frame exists at each point of a QC manifold. 
With respect to a qc-normal frame the above divergence reduces to 

V*(hA) = -{V ea A){he a ). 

Lemma 4.1. Suppose (M, rj, Q) is a quaternionic contact manifold with constant qc-scalar curvature. 
For any function h we have the following formulas 



V* (J2dh^ s )I s A s j = J2 ^dh(I s e a ^ s )A s (e a ) 

s=l s=l 
3 3 

V* (Y,dh(£ s )I s A) = ]T Wdh(I s e a ,^)A(e a ). 



8=1 8 = 1 

Proof. Using the identification of the 3-dimcnsional vector spaces spanned by {£1 , £2 , £3} and {I\ ,12,1$} 
with M 3 , the restriction of the action of Sp(n)Sp(l) to this spaces can be identified with the action 
of the group 50(3). With this in mind, one verifies easily that the 1-forms A, Y^ s =i dh(t; s )I s A s = 
— Yli=i dh(£i)[£j, £k] an d E s =i ^fc)-^^ are Sp(n)Sp(l) invariant on H. Thus, it is sufficient to 
compute their divergences in a qc-normal frame. To avoid the introduction of new variables, in this 
proof, we shall assume that {ei, . . . , e4„, £1, £2, £3} is a qc-normal frame. 

We apply (2.14). Using that the Biquard connection preserves the splitting of TM, we find 

V*[a,6] = -5(V eo (T(a,6)),e„) 

= - 5 ((V eo T)(a,6),e a ) - <?(T(V eo a,6),e ) - g(T(^, V e „6), **)■ 
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From Bianchi's identity we have (<ja.b,c means a cyclic sum over (A, B,C)) 
ff((V e .T)(&,6),e ) = -<?((V 6 T)(k,e ),e ) - fl ((V fo T)(e ,fi),e„) 

- .9(c ea ,( 3 \r{T(e a , £1), &)} , e a ) + g{a Sa , fl , ?2 |i?(e a , fi)£ 2 }, e a ) 

= -g(T(T(e a ,a),6),ea) - ff(T( T(£i, 6), e a ), e a ) - ff(T(T(6, e„), 6), e„) 

= 5(T(r(a,e a ),6),e ) - g(T(Tfa, e a ),^), e a ) - g{T(Tfa,&),e a ),e a ), 

taking into account that as mappings on H the torsion tensors T(£i,X) and the curvature tensor 
-ft(£i>£2) are traceless, so 5((V^ 1 T)(^2 ) e a ), e a ) and g(R(£i, $,%)e a , e a ) = 0, while the connection 
preserves the splitting, to obtain the next to last line. The last term is equal to zero as 

Seal Seal 
5 (T(T(£i,£ 2 ),e a ),e a ) = g(T{- — ^ - & - [fr, g 2 ] g , e a ), e a ) = - ^ - 2 ^ g(T(g 3 , e a ), e a ) = 0, 

taking into account that the torsion T^ 3 is traceless and T([£i, ^2]//, e a ) is a vertical vector. On the 
other hand, 

g(T{ T iiZa,£,2),e a ) - #(T(T£ 2 e a ,£i),e ) 

ff(r(e 6 ,£2),e a )ff(T(£i,e ),e6) - g(T(e 6 ,^),e ) g(T(£ 2 , e a ),e b ) 

= g(T(&,e b ),e a )g(T(£ 1 ,e a ),e b ) - g{T(£i, e b ), e a ) g(T(&, e a ), e 6 )J = 0. 

The equalities V*(/iA 1 ) = V*(/2^4 2 ) = with respect to a qc-normal frame can be obtained 
similarly. Hence, the first formula in Lemma 4.1 follows. 

We are left with proving the second divergence formula. Since the scalar curvature is constant 
(2.12) implies 

3 

(4.1) A(*) = -2 

s=l 

Fix an s € {1,2,3}. Working again in a qc-normal frame we have 

3 

(V ea A)(I s e a ) = -2^(V eo p t )(/ s e ai &). 
t=i 

A calculation involving the expressions (2.11) and the properties of the torsion shows that 

(4-2) tr(p t o/ s ) = - L— ktScaL 

2[n + I) 

The second Bianchi identity 

= .g((V ea i?)(/ s e a ,^)e b , I t e b ) + .g((V /sea i?)(£ t , e a )e b , I t e b ) + .g((V ?t i?)(e a , I s e a )e b , I t e b ) 

)e b ,I t e b ) + g(R(T(I s e a ,£t),e a )e b ,Ite b ) + g(R(T(£ t , e a ), I s e a )e b , I t e b ), 

together with the constancy of the qc-scalar curvature and (4.2) show that the third term on the 
right is zero and thus 
3 

^{2(V eoft )(7 s e a ,6) - 2p t (T(t t ,I s e a ),e a ) + p t (T(e a ,I s e a ),^)] = 0. 
t=i 

Substituting (2.5) in the above equality we come to the equation 

33 3 
(4.3) ^(V eaPt )(^e a ,6) = ^pt(T(6,/ s e a ),e o )-4n^p t (6,6)=0, 
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where the vanishing of the second term follows from (2.13), while the vanishing of the first term is 
seen as follows. Using the standard inner product on End(H) 

-in 

g(C, B) = tr(B*C) = £ g(C(e a ), B(e a )), 

a=l 

where C,B £ End(H), {ei, e4„} is a g-orthonormal basis of H , the definition of T® , the formulas 
in Theorem 2.5 and Proposition 2.3 imply 

3 

= aiPi^^h) + g(p2,T° 2 h) + ff(p 3 ,r| 3 7i) - g(pi,u) - g(p 2 ,I 3 u) + g(p 3 ,I 2 u) 
= g( Pl ,Tlh) + gfaT&h) + gfaTgh) 

+ 9 (2 { T^-J 1 - 2I 2 u - ^L_I 2 , T° 2 I l)+ 5 (2(T0)- + -/ 2 - 2/ 3 . - ^- h , 

= -2 5 ((T°)-+J 2 ,T°) + 2 5 ((T| 3 )+-,T°) + 2g((Tl)-+-I 3 ,Tl) 

= 2 5 ((T| 3 ) + - ) (T0)+-) + 2g((Tl)-+-,I 3 (Tl)+-) 

= 2g{I 2 {Tl)+-,I 2 (Tl)+-) 2g{h{Tl)- + -,I 2 {Tiy-) = 0. 

Renaming the almost complex structures shows that the same conclusion is true when we replace I\ 
with I 2 or I 3 in the above calculation. 

Finally, the second formula in Lemma 4.1 follows from (4.1) and (4.3). □ 

We shall also need the following one-forms 

L>iP0 = -hr 1 ^'' (X,Vh) 

(4.4) D 2 (X) = -/T 1 T (r+ ~ (X,V/i) 

D 3 (X) = -/r 1 T°" + (X,V/i) 

For simplicity, using the musical isomorphism, we will denote with D 1; D 2 , D 3 the corresponding 
(horizontal) vector fields, for example g{D ll X) = Di(X) \/X G H. Finally, we set 

(4.5) D = L>! + D 2 + D 3 = -h^T^X^h). 

Lemma 4.2. Suppose {M,rf) is a quaternionic contact manifold with constant qc-scalar curvature 
Seal = 16n(n + 2). Suppose fj = jf^rj has vanishing [— l]-torsion component T° = 0. We have 
^3 3 
D(X) = -h- 2 (3Wdh(X,Wh) - Vdh(I s X,I s Vhf) + h- 2 ^2dh(£ s ) dh(I s X). 

s=l s=l 

and the divergence of D satisfies 

V*L> = |T°| 2 -hr 1 g{dh,D) - hr 1 (n + 2)g{dh 1 A). 

Proof, a) The formula for D follows immediately from (3.2). 

b) We work in a qc-normal frame. Since the scalar curvature is assumed to be constant we use 
(2.15) to find 

V*D = -h- 1 dh{e a )D{e a ) - h- x V*T°(Vh) - h~ x T Q {e a , e b ) Wdh(e a , e b ) 

= -h- 1 dh(e a )D(e a ) - fc-^n + 2) dh(e a )A(e a ) - g{T°, h' 1 Vdh) 

= |T°| 2 - h- x dh{e a )D{e a ) - hT^n + 2) dh(e a )A(e a ), 
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using (3.2) in the last equality. □ 
Let us also consider the following one-forms (and corresponding vectors) 
F S (X) = -h^J^^XJsVh), XeH s = 1,2,3. 
From the definition of F\ and (4.4) we find 

Ft(X) = -h- l T°(X,IiVh) 

= -h^T^'iXJxVh) - h- 1 T^ + '(X,I 1 Vh) - h,- 1 T°~~ + (X, hVh) 

= h' 1 ^" (hX,Wh) - h^T ^' (hX,Vh) - h^T ''* (hX,Vh) 

= -Di(hX) + D 2 {hX) + D 3 (hX). 
Thus, the forms F s can be expressed by the forms D s as follows 

Fi{X) = -DtihX) + D 2 {hX) + D 3 (hX) 

(4.6) F 2 (X) = DxihX) - D 2 {I 2 X) + D 3 (I 2 X) 

F 3 {X) = D x {hX) + D 2 (I 3 X) - D 3 (I 3 X). 

Lemma 4.3. Suppose (M, 77) is a quaternionic contact manifold with constant qc-scalar curvature 
Scal = 16n(n + 2). Suppose fj = j^rj has vanishing [— l]-torsion component, T = 0. We have 

3 3 

V* (j2 dh (Zs)F s ) = Vdh(I s e a ,Z a )F a (I a e a ) 

s=l s=l 

3 

+ h- 1 Y\dKis)dh{I s e a )D{e a ) + (n + 2) dh(£ s )dh(I s e a ) A(e a ) . 

s=l 

Proof. We note that the vector X^ s =i dh(£ s )F s is an Sp(n)Sp(l) invariant vector, hence, we may 
assume that {ei, . . . , e± n , £i, £21 £3} is a qc-normal frame. Since the scalar curvature is assumed to 
be constant we can apply Theorem 2.6, thus V*T° = (n + 2)A Turning to the divergence, we 
compute 

33 3 

(4.7) ^(^dh^F^ = J2[vdh(e a ,UFs(ea)] - ^hT 1 dh(Z a )VT°(I 8 Vh) 

s=l s=l s=l 

3 

+ Y,[h~ 2 dh(£ s )dh(e a )T°(e a ,I s e b )dh{e b ) - h' 1 dh(£ s )T°(e a , I s e b ) X7dh(e a ,e b ) 

s=l 

3 3 

= J2[ Vdh(e a ,Z s )F s (e a )\ - Y, h ~ ldh (ts)V*T°(I s Vh) 

s=l s=l 
3 

+ h- 1 dh{t- a ) dh(I s e a ) D{e a ) 

s=l 

3 

= Y\ Vdh(e a ,£ s )F a (e a ) + h^ 1 dh(£ s ) dh(I s e a ) D(e a ) + hT^n + 2) dhfa) dh(I s e a ) A(e a ) 

s=l 

using the symmetry of T° in the next to last equality, and the fact 

3 

T°(e a ,he b ) Vdh(e a ,e b ) = - /i _1 V 'd/i[ SJ/m ][-i](e , he b ) Vdh\ sym ] (e , e b )~y^ dh(£ s )u s (e a , e b ) 



s=l 



0, 



which is a consequence of (3.2), the formula for the symmetric part of \7dh given after (3.3) and the 
zero traces of the [-l]-componcnt. Switching to the basis {I s e a : a = 1, . . . ,4n} in the first term of 
the right-hand-side of (4.7) completes the proof. □ 
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At this point we restrict our considerations to the 7-dimensional case, i.e. n = 1. Following is 
our main technical result. As mentioned in the introduction, we were motivated to seek a divergence 
formula of this type based on the Riemannian and CR cases of the considered problem. The main 
difficulty was to find a suitable vector field with non-negative divergence containing the norm of 
the torsion. The fulfilment of this task was facilitated by the results of [IMV], which in particular 
showed that similarly to the CR case, but unlike the Riemannian case, we were not able to achieve 
a proof based purely on the Bianchi identities, see [IMV, Theorem 4.8]. 

Theorem 4.4. Suppose (A/ 7 , 77) is a quaternionic contact structure conformal to a 3-Sasakian 
structure (M r ,fj), V = Jh V- If Seal v = Scal 7 ~ t = 16n(n + 2), then with f given by 

f =\ + h + \h~ l \Vh\ 2 , 

the following identity holds 

V*(/£> + ^2dhfo)F a + 4j2dh({ s )I s A 3 - ^-J2 dh &) I * A )= f\ T °\ 2 + h(QV,V). 

s=l s=l s=l 

Here, Q is a positive semi-definite matrix and V = (D\, D 2 , D$, A\, A 2 , A3) with A s , D s defined, 
correspondingly, in (3.5) and (4.4). 

Proof. Using the formulas for the divergences of D, X) s =i dh(£ s ) F s , J2 s =i dh(£ s )I s A s and Y^=i dh(£ s ) I S A 
given correspondingly in Lemmas 4.2, 4.3 and 4.1 we have the identity (n = 1 here) 

(4.8) V*(/£> + Y,dh(£ s )F s + Aj2dh(C s )I s A s - y £ d/»(&) J.a) 

S—1 S—1 S—1 

= (dh(e a ) - h- 2 dh(e a )\Vh\ 2 + hi- 1 Vdh(e a ,Vh)^j D(e a ) 

+ f(\T°\ 2 -h- l dh{e a )D{e a ) - fc-^n + 2) dh{e a )A(e a )) 

3 3 
+ Vdh(I s e a ,Z s )F s (I s e a ) + h^^dh^) dh(I s e a ) D(e a ) + (n + 2) dh(£ a ) dh(I a e a ) A(e a ) 



S = l 8=1 



3 10 3 

+ lJ2Vdh(I s e a ,£ s )A s (e a ) -- ^ Vdh(I s e a ,£ s ) A(e a ) 

s=l s=l 

1 1 3 

= (dh(e a ) - -h- 2 dh{e a )\Vh\ 2 + -h- 1 Vdh(e a ,Vh)) ^A(e«) 

t=i 

3 3 
+ /(|T°| 2 -h- l dh{e a ])(^D t {e a )) - fh- 1 {n + 2)dh{e a )(^A t {e a )) 

t=i t=i 
+ Vdh(I 1 e a ,^ 1 ){D 1 (e a ) - D 2 (e a ) - D 3 (e a j) + Vdh(I 2 e a , &) (--Di(e a ) + D 2 (e a ) - D 3 (e a )) 

+ Vd/i(/ 3 ea,6)(-^i(ea) - D 2 (e a ) + D 3 (e a )) 
33 33 
+ h- x (Y,dh{t s )dh{I a e a )) (X)A(ea)) + h- 1 (n + 2)(Y / dh(£ s )dh(I s e a ))(Y / A t (e a )) 

s=l t=l s=l t=l 

3 10 3 3 

+ 4^Vd/i(/ s e ,^)A 8 (e ) -y£ Vd/*(I.e a ,6)) E ^t(e„)), 

s=l s=l t=l 

where the last equality uses (4.6) to express the vectors F s by D s , and the expansions of the vectors 
A and D according to (3.5) and (4.5). Since the dimension of M is seven it follows U — U = 
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[ Vdh — 2h 1 dh g) d/i]™™ = 0. This, together with the Yamabe equation (3.4) , which when n = 1 
becomes Ah = 2 - Ah + S^V/il 2 , yield the formula, cf. (3.13), 

3 

(4.9) Vdh{X,Vh) + ^2Vdh(I s X,I s Vh) - ( 2 - Ah + Zh- l \Vh\ 2 ) dh{X) = 0. 

s=l 

From equations (4.4) and (3.12) we have 

D X {X) = h~ % dh{ix) dh(hX) + -h~ 2 [Vdh{X,Vh) + Vdh(hX,hVh) 

4 

- Vdh{I 2 X,I 2 Vh) - Wdh(I 3 X,I 3 Wh)], 
D 2 {X) = h- 2 dh(&) dh(I 2 X) + -^hr 2 [S7dh (X, V/i) - \Jdh{hX,hVh) 

+ Vdh(I 2 X,I 2 Vh) - Vdh(I 3 X,I 3 Vh)}, 
D 3 (X) = h- 2 dh{£ 3 ) dh{I 3 X) + ^h~ 2 [\7dh(X 7 \7h)~ Vdh(hX,hVh) 

- Vdh(I 2 X,I 2 Vh) + Vdh(I 3 XJ 3 Vh)]. 

Expressing the first term in (4.9) by the rest and substituting with the result in the above equations 
we come to 

(4.10) Di{e a ) = hi- 2 {2 - Ah + Sh'^hf) dh(e a ) + ^ 2 dhfe) dh(I t e a )+ 

X -h- 2 [- Vdh(Ijea,IjVh) -Vdh{I k e a ,I k Vh)] . 

At this point, by a purely algebraic calculation, using Lemma 3.1 and (4.10) we find: 

22 , 2 , 2 , 11 1 1 

Y A, - -A 2 - -A 3 + y Dx - -D 2 - -D 3 



-3h-i [ 1 2 



h-'dhiea) + h- 2 \Vh\ 2 ^j dh(e a ) + 3h- 2 dh (Cs) dh{I s e a )^j 

+ Ih-^dhihe*,^) - ^h- l Vdh(I 2 e a ,&) - ^h- l Vdh{he a ,&). 

Similarly, 

3A X - A 2 - A 3 + 2Di = (-2k- 1 + \^r 2 + /^ 3 |V/i| 2 ^ dh(e a ) - ^hr 2 Vdh(I s e a , I s Vh) 

3 

+ Vdfc(iie ,fc) - h- 1 Vdh(I 2 e a ,&) - h- 1 Vdh(I 3 e a ,Z 3 ) + h- 2 Y,dh(£ s )dh(I s e a ). 
On the other hand, the coefficient of Ai(e a ) in (4.8) is found to be, after setting n = 1, 

-sfl + \ h ~ 1 + ^h- 2 \S7h\Ah- 1 dh(e a )+ 3h- 2 dh&) dh(I s e a )) 



Zh- 1 Vdh(I 1 e a ,£ 1 ) - Y h ~ l V dh V2e a ,&) - y^ _1 V^(/ 3 ea,6) 
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while the coefficient of Di(e a ) in (4.8) is 
(4.11) dh(e a ) 



jh~ 2 dh(e a )\Vh\ 2 + X -hT x Vdh(e a ,Vh)- fh- l dh{e a ) 



Wdh(he a ,^ 



Vdh(I 2 e a ,&) - V^(/ 3 e a ,e 3 )^i(ea) 



h-^Y^dh^dhilsea)). 



Substituting Vdh(e a ,Vh) according to (4.9), i.e., Vd/i(e a ,V/i) = — Y^s=i ^dh (I s e a , I s Vh) + 
(2-4/i + 3h~ 1 \S7h\ 2 ) dh(e a ) and using the definition of / transforms the above expression into 



dh(e a ) - hi- 2 dh{e a )\Vh\ 2 - Q + h + ih _1 |Vft| 2 j h^dJi^a) 

+ hi- 1 (- ^~yVdh(I s e a ,I s Vh) + ( 2 - Ah + Sh^lVhl 2 ) dh(e a )j 

+ Vdh(/ie ,fi) - Vdfc(J a e ,6)- Vdfc(J 3 e.,6)i3i(e a ) + ^ ( ^ d/i(&) d/i(/ s e ) 

\ s =i y 

Simplifying the above expression shows that the coefficient of D\{e a ) in (4.8) is 
-2 + hi- 1 + /i" 2 |V/i| 2 ) dfc(e„) - \h~ x ^ Vd/i (I s e a ,I s Vh) ^ 

+ Vdh(he a ,^) - Vd/i(/ a e a ,£ a )- Vdh(I 3 e a: &) + fr- 1 ( ]T d/i(&) d/i(/ s e a ) 



s=l 



Hence, we proved that the coefficient of Di(e a ) in (4.8) is /i (3 A\ — — A 3 + 2Di) (e a ), while 
those of A 1 (e a ) is A\ — § j4 a — § A3 + ^ D\ — \ D2 — § Z? 3 )(e a ). A cyclic permutation 
gives the rest of the coefficients in (4.8). With this, the divergence (4.8) can be written in the form 



V*(/D + J2 dh tis)Fs + Aj2dh($ s )I s A s - ^-^dh^.) I,A 



f\T° 



h<T h2 . 3 {g(D 1 ,3A 1 - A 2 - A 3 + 2D, 



, 22 
g(Ai, —Ai 



2 , 2 A 11 1 1 

- A 2 - - A 3 + — Di - -D 2 - -D 3 

3 3 3 3 3 



where (71,2.3 denotes the sum over all positive permutations of (1,2,3). Let Q be equal to 



Q~ 
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so that 

3 in 3 
V*(/D + ]Td/i(&)F a + ^dh{^)I s A s - —Y^M^IsA) = f\T°\ 2 + h(QV,V), 

s=l s=l 

with V = {Di, D2, -D3, Ai, A2, A3). It is not hard to see that the eigenvalues of Q are given by 

{0, 0, 2(2 + V2), 2(2-V2), 10, 10}, 
which shows that Q is a non-negative matrix. □ 

5. Proofs of the main theorems 

The proofs rely on Theorem 4.4 and the following characterization of all qc-Einstcin structures 
conformal to the standard qc structures on the Heisenberg group. 

Theorem 5.1. [IMV, Theorem 1.2] Let = ^-0 be a conformal deformation of the standard Re- 
structure on the quaternionic Heisenberg group G (H). Then is qc-Einstein if and only if, up 
to a left translation the function h is given by 

(5.1) h = c [(1 + y\q\ 2 f + v 2 (x 2 + y 2 + z 2 )] , 

where c and v are any positive constants. 

Consider first the case of the (seven dimensional) sphere. 

5.1. Proof of Theorem 1.1. Integrating the divergence formula of Theorem 4.4 we see that ac- 
cording to the divergence theorem established in [ M\ , Proposition 8.1] the integral of the left-hand 
side is zero. Thus, the right-hand side vanishes as well, which shows that the quaternionic contact 
structure rj has vanishing torsion, i.e., it is also qc-Einstcin according to Theorem 2.4. 

Next we bring into consideration the 7-dimensional quaternionic Heisenberg group and the quater- 
nionic Cayley transform as described in [IMV, Section 5.2]. The quaternionic Heisenberg group of 
dimension 7 is G (H) = H x ImH. The group law is given by (q',oj') = (q ,u> ) o (q,ui) = 
(q + q, uj + uj + 2 Im q q), where q, q £ EI and lo,lo G ImH. The left-invariant orthonormal 
basis of the horizontal space is 

Ti = — +2x 1 — + 2y 1 — + 2z 1 —, X 1 = 2t 1 — - 2z Y — + 2y 1 — 

oti ox oy oz 0x1 ox oy Oz 

d „ , 3 „, 3 „ , 3 „ d „ , 8 „ , 3 .3 

Y i = a- + 2 * a"" 2 * a- _2a! V' Zl = a 2 ^V + 2x V - 2 * V 

oyi ox oy oz oz\ ox oy Oz 

using q = t\ + ix\ + j yi + kz\ and lu = ix + jy + kz. The central (vertical) orthonormal 
vector fields £1,^2, £3 are described as follows 

d d d 

6 = 2— 6 = 2^ 6 = 2^-. 

ax oy Oz 

Let us identify the (seven dimensional) group G (H) with the boundary E of a Siegel domain in 

H x H, 

E = {(<?', p') eHxH : Kp' = | g ' 2 }. 

E carries a natural group structure and the map (q,u>) i— > (g, |g| 2 — w) G E is an isomorphism 
between G (H) and E. 

The standard contact form, written as a purely imaginary quaternion valued form, on G (M) 
is given by 20 = (cko — q ■ dq + dq ■ q), where • denotes the quaternion multiplication. 
Since dp = q ■ dq + dq ■ q — du>, under the identification of G (H) with E we also have 
20 = —dp' + 2dq' ■ q'. Taking into account that is purely imaginary, the last equation can be 
written also in the following form 

4 = (dp' - dp') + 2dq'-q' - 2q'-dq'. 
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The (quaternionic) Cayley transform is the map G : S \ {(—1,0)} i— > E from the sphere S = 
{(q,p) eHxI : |g| 2 + \p\ 2 = 1} c H x H minus a point to the Heisenberg group E = {(qi,pi) £ 
H x H : Upi = |9i| 2 }, with 6 defined by 

(5.2) ( qi , Pl ) = G ((q,p)), q^il+py'q, Pl = (1+p)- 1 (1-p). 
with an inverse (q, p) = C _1 f(5i,£>i)J given by 

(5.3) q = 2(l+p 1 )- 1 q u p = (l- Pl )(l +Pl )-\ 

The Cayley transform is a conformal quaternionic contact diffcomorphism between the quaternionic 
Heisenberg group with its standard quaternionic contact structure O and S \ {(—1,0)} with its 
standard structure fj, see [IMV], 

(5.4) a ■ a = ir^e, 

where A = | is a unit quaternion and 77 is the standard quaternionic contact form on the 
sphere, 77 = dq ■ q + dp ■ p — q ■ dq— p ■ dp. Hence, up to a constant multiplicative factor and 
a quaternionic contact automorphism the forms 6*77 and O are conformal to each other. It follows 
that the same is true for 6*77 and 0. In addition, is qc-Einstein by definition, while 77 and hence 
also 6*77 are qc-Einstein as we observed at the beginning of the proof. According to Theorem 5.1, up 
to a multiplicative constant factor, the forms 6*77 and 6*77 are related by a translation or dilation on 
the Heisenberg group. Hence, we conclude that up to a multiplicative constant, 77 is obtained from 
77 by a conformal quaternionic contact automorpism which proves the first claim of Theorem 1.1. 
From the conformal properties of the Cayley transform and [Val, Va] it follows that the minimum 
X(S 4n+3 ) is achieved by a smooth 3-contact form, which due to the Yamabe equation is of constant 
qc-scalar curvature. This shows the second claim of Theorem 1.1. 



5.2. Proof of Theorem 1.3. Let D 1,2 be the space of functions u £ L 2 (G (H)) having distribu- 
tional horizontal gradient |Vw| 2 = \Tiu\ 2 + \Xiu\ 2 + \Yiu\ 2 + \Z\u\ 2, £ L 2 (G (H)) with respect to the 
Lebesgue measure dH on R 7 , which is the Haar measure on the group. Let us define the constant 
(2* = 5/2 here) 




Let v be a function for which the infimum is achieved. Note that such function exists by [Val] or 
[Va]. Furthermore, A = S^ 2 , where S2 is the best constant in the L 2 Folland-Stcin inequality (1.1), 
since v £ D 1 - 2 implies |u| £ D 1,2 and the gradient is the same a.e.. From the choice of v we have 

A = J |Vv| 2 dH, J v r dH = 1. 

G (H) G (H) 

Writing the Euler-Lagrangc equation of the constrained problem we see that v is a non-negative 
entire solution of (T- 2 + X\ + Y? + Zf)v = -A v 3 / 2 . By [GV2, Lemma 10.2] ( see [Va] or [Val, 
Theorem 10.3] for further details) v is a bounded function. Similarly to [FSt, Theorem 16.7] it follows 
v is a Lipschits continuous function in the sense of non- isotropic Lipschits spaces [F]. Iterating this 
argument and using [F, Theorem 5.25] we see that v is a C°° smooth function on the set where it is 
positive, while being of class T^, the non-isotropic Lipschits space, for some p > 0. In particular 
v is continuously diffcrcntiable function by [F, Theorem 5.25]. Applying the Hopf lemma [GV1, 
Theorem 2.13] on the set where v is positive shows that v cannot vanish, i.e., it is a positive entire 
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solution to the Yamabe equation. The positivity can also be seen by the Harnack inequality, see [W] 
for example. Let u = A 2 *- 2 v, then u is a positive entire solution of the Yamabe equation 
(5.5) (T 2 + X\ + Y 2 + Z\)u = -u 3/2 

^From the definition of u, we have 

A = ( J \Vu\ 2 dff)* = ( J u 5 ' 2 dHy. 

G (H) G (H) 

We shall compute the last integral by determining u with the help of the divergence formula. 

As before, let be the standard contact form on G (H) identified with E. Using the inversion 
and the Kelvin transform on G (H), cf. [GV2, Sections 8 and 9], we can see that if = t^-0 has 
constant scalar curvature, then the Cayley transform lifts the qc structure defined by to a qc 
structure of constant qc-scalar curvature on the sphere, which is conformal to the standard. The 
details are as follows. Let us define two contact forms 0i and 02 on E setting 

0! = U 4 /W- 2 >0, and 2 = (Xu) 4 /^-V £- *L, 

\P\ \P\ 

where u is as in (5.5), Xu is its Kelvin transform, see (5.8) for the exact formula, and Q is the 
homogeneous dimension of the group. Notice that At t^j defines the same qc structure on the 

Heisenbcrg group as and Xu is a smooth function on the whole group according to [GV2, Theorem 
9.2]. We are going to see that using the Cayley transform these two contact forms define a contact 
form on the sphere, which is conformal to the standard and has constant qc-scalar curvature. 

Let Pi = ( — 1,0) and P 2 = (1,0) be correspondingly the 'south' and 'north' poles of the unit 
sphere S = {\q\ 2 + \p\ 2 = 1}. Let Ci and 62 be the corresponding Cayley transforms defined, 
respectively, on S\ {Pi} and S \ {^2}- Note that Ci was defined in (5.2), while C2 is given by 



(5.6) (q 2lP2 ) = C 2 (iq,p)), q 2 = -(1 - p)' 1 q, p 2 = (1 - pV (1 + p). 
In order that 0i and ©2 define a contact form 77 on the sphere it is enough to see that 

(5.7) e^.g) = 2 oe 2 oei- 1 ( P , (? ), i.e., ©! = (e 2 oe 1 - 1 )*e 2 . 

A calculation shows that C 2 o Cj" 1 : E — > E is given by 



Q2 = -p 1 1 qi, P2 = Pi \ 

or, equivalently, in the model G (H) 

<\ |2 \-l w i 

92 = -{ml -Wl) qi, U 2 = - 1 14, I — 12- 

m\ + Fir 

Hence, a = C 2 o C^ 1 is an involution on the group. Furthermore, with the help of (5.4) we calculate 

1 -r, Pi 



61,062 = 1 pr^lQfi, fJi 

\Pl\ \P 



which proves the identity (5.7). Using the properties of the Kelvin transform, [GV2, Sections 8 and 
9], 

(5.8) (Xu)(q',p') = f \p'-( Q - 2)/2 u(a( q >,p')), 

we see that u and %u are solutions of the Yamabe equation (5.5). This implies that the contact 
form j] has constant qc-scalar curvature, equal to 4 q^" 2 ^ • 

Notice that r\ is conformal to the standard form 77 and the arguments in the preceding proof imply 
then that r\ is qc-Einstein. A small calculation shows that this is equivalent to the fact that if we 
set 
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(5.9) u = 2 10 [(1 + M 2 ) 2 + |c| 2 ]- 2 , 

then u satisfies the Yamabe equation (5.5) and all other nonnegative solutions of (5.5) in the space 
D 1 ' 2 arc obtained from u by translations and dilations, 

def _ 

(5.10) i~(q otUJo )U (q, uj) = u(q Q + q,u + w ), 

(5.11) u\(q) d = A 4 u(Ag,A 2 w), A > 0. 

Thus, u which was defined in the beginning of the proof is given by equation (5.9) up to transla- 
tions and dilations. This allows the calculation of the best constant in the Folland-Stein inequality, 
see [GV1, (4.52)], 

A 5 = rdH = 2 25 tt 7 / 2 - W 



[(i + M 2 ) 2 + M 2 )] 5 r(7) 12 

G (H) 

where T is the Gamma function. Hence 

c _ a -1/2 _ 2 ^ 

- ^375- 

Recalling the relation between u and v we find that the extremals in the Folland-Stein embedding 
are given by 

» - ^[(i + l^l 2 ) 2 + M 2 r 2 

and its translations and dilations. The proof of Theorem 1.3 is complete. 
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